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November 7, 2024 at 12:27 


1. Intro. This is a modification of the program SSXCC, extending it by adding proof logging. We will 
use the VeriPB program as the external verifier. This program creates two extra output files: an OPB file 
containing a system of linear inequalities to model our XCC problem, and a PBP file which contains our 
proof logging. The proof logging extensions were introduced by Filip Stappers in 2024. 

I shall accept the DLX input format used in the previous solvers, without change, so that a fair comparison 
can be made. (See the program DLX2 for definitions. Much of the code from that program is used to parse 
the input for this one.) 

My original attempt, SSXCO, kept the basic structure of DLX1 and changed only the data structure 
link conventions. The present version incorporates new ideas from Christine Solnon’s program XCC-WITH- 
DANCING-CELLS, which she wrote in October 2020. In particular, she proposed saving all the active set 
sizes on a stack; program SSXCO recomputed them by undoing the forward calculations in reverse. She also 
showed how to unify “purification” with “covering.” 

In August, 2023, Christine told me about two further improvements: We can easily recognize most cases 
where an item has only one option left (a “forced move”). And in such cases, it isn’t necessary to save the 
domain sizes of the other active elements, because we won’t need that information when backtracking. 


2. After this program finds all solutions, it normally prints their total number on stderr, together with 
statistics about how many nodes were in the search tree, and how many “updates” were made. The running 
time in “mems” is also reported, together with the approximate number of bytes needed for data storage. 
(An “update” is the removal of an option from its item list, or the removal of a satisfied color constraint 
from its option. One “mem” essentially means a memory access to a 64-bit word. The reported totals don’t 
include the time or space needed to parse the input or to format the output.) 


7##define o memst+ /* count one mem */ 

##define 00 mems += 2 /* count two mems */ 

7##define 000 mems += 3 /* count three mems */ 

#define O "%" /* used for percent signs in format strings «/ 

#define mod % /* used for percent signs denoting remainder in C */ 

##define maz_level 5000 /* at most this many options in a solution */ 

#define maz_cols 100000 /* at most this many items */ 

7##define maz_nodes 10000000 /* at most this many nonzero elements in the matrix */ 
#define savesize 10000000 /*x at most this many entries on savestack */ 


#tdefine bufsize (9 * maz_cols + 3) /* a buffer big enough to hold all item names */ 
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Here is the overall structure: 


#include <stdio.h> 
#include <stdlib.h> 
#include <string.h> 
#include <ctype.h> 
#include "gb_flip.h" 


typedef unsigned int uint; /* a convenient abbreviation */ 

typedef unsigned long long ullng; /* ditto «/ 

(Type definitions 8 ); 

(Global variables 4); 

(Subroutines 11); 

main(int argc,char xargv[]){ register int c, cc, i, j, k, p, pp, qg, r, 8, t, cur-chotce, cur_node, 


best_itm, line , optnr; 

(Process the command line 5); 
(Input the item names 15); 
(Input the options 17); 
(Output our problem to the OPB file for proof logging 27); 
(Output header of the PBP file for proof logging 32); 
if (vbose & show_basics) (Report the successful completion of the input phase 24); 
if (vbose & show_tots) (Report the item totals 25); 
imems = mems,mems = 0; 
if (baditem) (Report an uncoverable item 23 ) 
else { 

if (randomizing) {Randomize the item list 26); 

(Solve the problem 38 ); 


} 


done: (Write contradiction line for proof logging 36 ); 


if (vbose & show_profile) (Print the profile 55 ); 
if (vbose & show_-max_deg ) 
fprintf (stderr , "The, maximum branching, degree, was,"O"d.\n", maxdeg); 
if (vbose & show_basics) { 
forintf (stderr, "Altogether,"O"1llu,solution"O"s,"O"llut+"O"llu,mems,", count, 
count =1?"":"s",imems,mems); 
bytes = (itemlength + setlength) * sizeof (int) + last_node * sizeof 
(node) + 2 * maal « sizeof (int) + maxsaveptr * sizeof (twoints); 
fprintf (stderr, "\"O"1llu updates, "O"1llu,bytes,"O"llunodes,", updates, bytes, nodes); 
fprintf (stderr, "\Wccosty"O"11d%%.\n", mems ? (200 * cmems + mems)/(2 * mems) : 0); 


me 


f (sanity_checking) fprintf (stderr, "sanity_checking,was,on!\n"); 
Close the files 6); 


—— 
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4. You can control the amount of output, as well as certain properties of the algorithm, by specifying 
options on the command line: 


e ‘v(integer )’ enables or disables various kinds of verbose output on stderr, given by binary codes such as 
show-choices ; 

e ‘m( integer )’ causes every mth solution to be output (the default is mO, which merely counts them); 

e ‘s(integer )’ causes the algorithm to randomize the initial list of items (thus providing some variety, 
although the solutions are by no means uniformly random); 

e ‘d(integer )’ sets delta, which causes periodic state reports on stderr after the algorithm has performed 
approximately delta mems since the previous report (default 10000000000); 

e ‘c( positive integer )’ limits the levels on which choices are shown during verbose tracing; 

e ‘C( positive integer )’ limits the levels on which choices are shown in the periodic state reports; 

‘1(nonnegative integer )’ gives a lower limit, relative to the maximum level so far achieved, to the levels 

on which choices are shown during verbose tracing; 

‘t (positive integer )’ causes the program to stop after this many solutions have been found; 

‘T( integer )’ sets timeout (which causes abrupt termination if mems > timeout at the beginning of a level); 

‘S( filename )’ to output a “shape file” that encodes the search tree. 

‘P( filename )’ (required) name of the “.opb file” and “.pbp file” for proof logging; 


7##define show_basics 1 /* vubose code for basic stats; this is the default */ 

#define show_choices 2 /* vbose code for backtrack logging */ 

##define show_details 4 /* vubose code for further commentary */ 

#define show_profile 128 /* ubose code to show the search tree profile «/ 

#define show_fullstate 256 /* vbose code for complete state reports */ 

7##define show_tots 512 /* vubose code for reporting item totals at start */ 

#define show_warnings 1024 /* vubose code for reporting options without primaries */ 
#define show_mazdeg 2048 /* vbose code for reporting maximum branching degree */ 


(Global variables 4) = 
int random_seed = 0; /* seed for the random words of gb_rand */ 
int randomizing; /* has ‘s’ been specified? */ 
int vbose = show_basics + show_warnings; /* level of verbosity */ 
int spacing; /* solution k is output if k is a multiple of spacing */ 
int show_choices.maz = 1000000; /* above this level, show_choices is ignored */ 
int show_choices.gap = 1000000; /* below level mazl — show-choices.gap, show_details is ignored */ 
int show_levels.max = 1000000; /* above this level, state reports stop */ 
int mazl; /* maximum level actually reached */ 
int mazsaveptr; /* maximum size of savestack */ 
char buf [bufsize]; /* input buffer «/ 
ullng count; /* solutions found so far */ 
ullng options; /* options seen so far */ 
ullng imems, mems, tmems, cmems; /* mem counts */ 
ullng updates; /* update counts «/ 
ullng bytes; /* memory used by main data structures */ 
ullng nodes; /* total number of branch nodes initiated */ 
ullng thresh = 10000000000; /* report when mems exceeds this, if delta £0 */ 
ullng delta = 10000000000; /* report every delta or somems */ 
ullng maxcount = *ffffffffffff{ffff; /* stop after finding this many solutions */ 
ullng temeout = *1ffffffff{fffffff; /* give up after this many mems */ 
FILE xshape_file; /* file for optional output of search tree shape */ 


char xshape_name; /* its name */ 
FILE xopb_file; /* file for OPB model of the problem we are solving */ 
char opb_name [255]; /* its name x/ 


FILE «pbp_file; /* file for our proof log */ 
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char pbp_name [255]; /* its name x/ 
int maxdeg; /* the largest branching degree seen so far */ 
int plLconstraints; /*x Current number of proof logging constraints written */ 


See also sections 9 and 39. 


This code is used in section 3. 


5. 


If an option appears more than once on the command line, the first appearance takes precedence. 


(Process the command line 5) = 
for (j = argc —1,k =0; 7; j——) 


switch (argv[j][0]) { 


case ’v’: k |= (sscanf (arguv[j] + 1,""O"da", &ubose) — 1); break; 

case ’m’: k |= (sscanf (arguv[j] +1,""O"d", &spacing ) — 1); break; 

case ’s’: k |= (sscanf (arguv[j] +.1,""O"d", &random_seed) — 1), randomizing = 1; break; 
case ’d’: k |= (sscanf (arguv[j] +1,""O"11d", &delta) — 1), thresh = delta; break; 

case ’c’: k |= (sscanf (argv|j] + 1,""O"d", &show_choices_maz ) — 1); break; 

case ’C’: k |= (sscanf (argv|j] + 1,""O"d", &show_levels.maxz) — 1); break; 

case ?1’: k |= (sscanf (argu[j] + 1,""O"d", &show_choices_gap ) — 1); break; 

case ’t’: k |= (sscanf (argv[j] + 1,""O"1l1d", &maxcount) — 1); break; 

case ’T’: k |= (sscanf (argv|j] + 1,""O"1lda", &timeout) — 1); break; 

case ’S’: shape._name = argu[j] + 1, shape_file = fopen (shape_name, "w"); 


if (ashape_file ) 
forintf (stderr, "Sorry, ,Ijcan’ tyopen file,‘"O"s’ for writing! \n", shape_name); 

break; 

case ’P’: strcpy(opb_name, argv[j] + 1), strcat(opb_name,".opb"), opb_file = fopen(opb_name, "w"); 
if (sopb_file) fprintf (stderr, "Sorry, ,Ican’ t open file, ‘"O"s’ for writing! \n", opbname); 
strepy (pbp_name, argv |j] + 1), strcat (pbp_name,".pbp"), pbp_file = fopen(pbp_name, "w"); 
if (=pbp_file) fprintf (stderr, "Sorry, ,Ican’ tyopen file, ‘"O"s’ for writing! \n", pbp_name); 
break; 

default: k = 1; /* unrecognized command-line option */ 


} 


if (k V sopb_file V apbp_file) { 


} 


fprintf (stderr, "Usage: ,"O"s,P<bar>, [v<n>] |, [m<n>] ,, [s<n>] ,[d<n>] ""[c<n>] ,[C<n>] , [1<n\ 
>] U[t<n>] ,[T<n>] , [S<bar>] |,<foo.d1lx\n", argu (0]); 
exit (—1); 


if (randomizing) gb_init_rand (random_seed ); 


This code is used in section 3. 


6. 
i 


f 
if 


(Close the files 6) = 
(shape_file) fclose(shape_file ); 
(opb_file) fclose (opb_file ); 


if (pbp_file) fclose (pbp_file); 


This code is used in section 3. 
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7. Data structures. Sparse-set data structures were introduced by Preston Briggs and Linda Torczon 
[ACM Letters on Programming Languages and Systems 2 (1993), 59-69], who realized that exercise 2.12 
in Aho, Hopcroft, and Ullman’s classic text The Design and Analysis of Computer Algorithms (Addison— 
Wesley, 1974) was much more than just a slick trick to avoid initializing an array. (Indeed, TAOCP exercise 
2.2.6—24 calls it the “sparse array trick.” ) 

The basic idea is amazingly simple, when specialized to the situations that we need to deal with: We can 
represent a subset S of the universe U = {2,%1,...,%n—1} by maintaining two n-element arrays p and gq, 
each of which is a permutation of {0,1,...,n — 1}, together with an integer s in the range 0 <s <n. In 
fact, p is the inverse of g; and s is the number of elements of $. The current value of the set S' is then simply 
{Lpo.+-+,Lp,_,}. (Notice that every s-element subset can be represented in s!(n — s)! ways.) 

It’s easy to test if x, € S, because that’s true if and only if gq, < s. It’s easy to insert a new element x, 
into S: Swap indices so that p, = k, q, = s, then increase s by 1. It’s easy to delete an element x; that 
belongs to S: Decrease s by 1, then swap indices so that p, = k and q, = s. And so on. 

Briggs and Torczon were interested in applications where s begins at zero and tends to remain small. 
In such cases, p and q need not be permutations: The values of ps, ps41, ---,; Pn—1 Can be garbage, and 
the values of gq, need be defined only when x, € S. (Such situations correspond to the treatment by Aho, 
Hopcroft, and Ullman, who started with an array full of garbage and used a sparse-set structure to remember 
the set of nongarbage cells.) Our applications are different: Each set begins equal to its intended universe, 
and gradually shrinks. In such cases, we might as well maintain inverse permutations. The basic operations 
go faster when we know in advance that we aren’t inserting an element that’s already present (nor deleting 
an element that isn’t). 

Many variations are possible. For example, p could be a permutation of {xo,21,...,2%n—1} instead of a 
permutation of {0,1,...,n —1}. The arrays that play the role of g in the following routines don’t have 
indices that are consecutive; they live inside of other structures. 
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8. This program has an array called item, with one entry for each item. The value of item[k] is an 
index x into a much larger array called set. The set of all options that involve the kth item appears in that 
array beginning at set[x]; and it continues for s consecutive entries, where s = size(a) is an abbreviation 
for set|a — 1]. If ttem[k] = x, we maintain the relation pos(x) = k, where pos(x) is an abbreviation for 
set[z — 2]. Thus item plays the role of array p, in a sparse-set data structure for the set of all currently 
active items; and pos plays the role of q. 

Suppose the Ath item x currently appears in s options. Those options are indices into nd, which is an 
array of “nodes.” Each node has three fields: itm, loc, and clr. Ifa <q < a+, let y = set|q]. This is 
essentially a pointer to a node, and we have nd[y].ttm = x, ndl[y].loc = q. In other words, the sequential 
list of s elements that begins at « = item|k] in the set array is the sparse-set representation of the currently 
active options that contain the kth item. The clr field nd[y].clr contains x’s color for this option. The itm 
and clr fields remain constant, once we’ve initialized everything, but the loc fields will change. 

The given options are stored sequentially in the nd array, with one node per item, separated by “spacer” 
nodes. If y is the spacer node following an option with ¢ items, we have nd[y].itm = —t. If y is the spacer 
node preceding an option with t items, we have nd[y].loc = t. 

This probably sounds confusing, until you can see some code. Meanwhile, let’s take note of the invariant 
relations that hold whenever k, g, x, and y have appropriate values: 


pos(item|k]) =k; nd[set[q]].loc =q; item[pos(x)] =x; set|[nd[y].loc] = y. 


(These are the analogs of the invariant relations p[q[k]] = q[p[k]] = & in the simple sparse-set scheme that 
we started with.) 

The set array contains also the item names. 

We count one mem for a simultaneous access to the itm and loc fields of a node. We don’t count mems 
for accesses to the linenr field, since these are only relevant for proof logging. 
#define size(x) set[(x) — 1] /* number of active options of the kth item, x */ 
#define pos(x) set[(ax) — 2] /* where that item is found in the item array */ 
#define Iname(x) set|(x) — 4] /* the first four bytes of x’s name */ 


#define rname(x)  set|(x) — 3] /* the last four bytes of x’s name */ 

#tdefine primextra 4 /* this many extra entries of set for each primary item */ 
7##define secondeztra 4 /* and this many for each secondary item */ 

#tdefine mazextra 4 /* maximum of primeztra and secondextra */ 


(Type definitions 8) = 
typedef struct node_struct { 


int itm; /* the item x corresponding to this node «/ 

int loc; /* where this node resides in x’s active set */ 

int clr; /* color associated with item x in this option, if any */ 

int linenr; /* line number from the input where this node comes from */ 
} node; 


See also section 10. 


This code is used in section 3. 
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9. (Global variables 4) += 


node nd[maz_nodes]; /* the master list of nodes */ 

int last_node; /* the first node in nd that’s not yet used */ 

int item|maz_cols]; /* the master list of items */ 

int second = maz_cols; /* boundary between primary and secondary items */ 
int last_itm; /* items seen so far during input, plus 1 */ 

int set[maz_nodes + maxextra * maz_cols]|; /* the sets of active options for active items */ 
int itemlength; /* number of elements used in item */ 

int setlength; /* number of elements used in set */ 

int active; /* current number of active items */ 

int oactive; /* value of active before swapping out current-choice items */ 
int baditem; /* an item with no options, plus 1 */ 

int osecond; /* setting of second just after initial input «/ 

int force|maz_cols]; /* stack of items known to have size 1 «/ 

int forced; /* the number of items on that stack «/ 


10. We're going to store string data (an item’s name) in the midst of the integer array set. So we’ve got 
to do some type coercion using low-level C-ness. 
( Type definitions 8) += 

typedef struct { 


int /, rT; 

} twoints; 

typedef union { 
unsigned char str[8]; /* eight one-byte characters */ 
twoints Ir; /* two four-byte integers */ 


} stringbuf; 
stringbuf namebuf ; 


11. (Subroutines 11) = 
void print_item_name (int k, FILE «stream) 


{ 


namebuf .lr.l = Iname(k), namebuf .lr.r = rname(k); 
forintf (stream, "\"O" .88", namebuf .str); 


See also sections 12, 13, 14, 28, 37, 44, 52, 53, and 54. 


This code is used in section 3. 
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12. An option is identified by the names of the items it contains. For proof logging we use the line number 
where the option appears in the input as the identifier of the option. Here is a routine that prints an option, 
given a pointer to any of its nodes. It also prints the position of the option in its item list. 
(Subroutines 11) += 
void print_option (int p, FILE stream) 
{ 
register int k, q, 2; 
x = nd[p].itm; 
if (p > last.node Va <0) { 
forintf (stderr, "Illegal,option,"O"d!\n", p); 
return; 
} 


for (q=p;; ) { 
print_item_name(a, stream); 
if (nd[q].clr) fprintf (stream, ":"O"c", nd[q].clr); 
qtr; 
x = nd[q].itm; 
if (1 <0) { 
q t= 2; 
fprintf (stream, "\,(Line,"O"d)", nd[q].linenr); 
x = nd[q].itm; 
} 
if (¢ =p) break; 
} 
k = nd[q].loc; 
fprintf (stream, "QC"O"dyofy"O"d) \n", k — x + 1, size(x)); 
} 
void prow(int p) 
{ 


print_option (p, stderr ); 


} 


13. When I’m debugging, I might want to look at one of the current item lists. 


(Subroutines 11) += 
void print_itm (int c) 
{ 
register int p; 
if (c < primextra V c > setlength V pos(c) < 0V pos(c) > itemlength V item[pos(c)| #c) { 
forintf (stderr, "Illegaljitem,"O"d!\n",c); 
return; 


forintf (stderr, "Item"); 
print.item_name(c, stderr); 
if (c < second) fprintf (stderr, "yC"O"dyofy"O"d) , ylength,"O"d:\n", pos(c) + 1, active, size(c)); 
else if (pos(c) > active) 
fprintf (stderr, "\,(secondary,"O"d, ypurified) , length,"O"d:\n", pos(c) + 1, size(c)); 
else fprintf(stderr, "\(secondary,"O"d) , ylength,"O"d:\n", pos(c) + 1, size(c)); 
for (p=c; p<c+ size(c); p++) prow(set[p)); 
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14. Speaking of debugging, here’s a routine to check if redundant parts of our data structure have gone 
awry. 
7##define sanity_checking 0 /* set this to 1 if you suspect a bug */ 


(Subroutines 11) += 
void sanity (void) 


register int k, x, 7, 1, r, q, 99; 
for (k =0; k < itemlength; k++) { 
x = item|k]; 
if (pos (i) #k) { 
forintf (stderr, "Bad. posyfield of item"); 
print.item_name (a, stderr); 


fprintf (stderr, "\("O"d,"O"d) !\n",k, x); 


for (¢ =0; i < last_node; i++) { 
1 = ndli].itm,r = nd{i].loc; 
if (1<0) { 
if (nd[it+r+1]).itm ¢ —r) fprintf(stderr, "Bad spacerin,nodes,"O"d, "O"d!\n",i,i+r+1); 
qq = 0; 
else { 
if (l>r) fprintf(stderr, "itm>loc in node,"O"d!\n", 4); 
else { 
if (set[r] £2) { 
fprintf (stderr, "Bad locyfield for option,,"O"d of yitem",r —1+ 1); 
print_item_name (1, stderr); 
fprintf (stderr, "\yinynode,"O"d! \n", 2); 


a 


if (pos(l) < active) { 
if (r <1+ size(l)) q=+1; else ¢=-1; /* in or out? */ 
if (q* qq <0) { 
forintf (stderr, "Flipped status_atoption,,"O"d of item",r —1+ 1); 
print_item_name (I, stderr); 
forintf (stderr, "\yinnode,,"O"d!\n", 2); 


qq = | 
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15. Inputting the matrix. Brute force is the rule in this part of the code, whose goal is to parse and 
store the input data and to check its validity. 

We use only four entries of set per item while reading the item-name line. 
#define panic(m) 

{ fprintf(stderr, ""O"s!\n"O"d:4"O".99s\n",m, p, buf); exit (—666); } 

(Input the item names 15) = 

line = 0; while (1) { 

if (fgets (buf , bufsize, stdin)) break; 

line ++ ; 
if (0, buf [p = strlen(buf) — 1] 4 ?\n’) panic("Inputylineway,jtoo, long"); 
for (p = 0; 0, isspace (buf [p]); p++) ; 
if (buf [p] =? |? V abuf [p]) continue; /* bypass comment or blank line «*/ 
last.itm = 1; 
break; } 
if (slast_itm) panic("No,items"); 
for (; 0, buf[p]; ) { 

o, namebuf .lr.l = namebuf .lr.r = 0; 

for (j =0; § <8A (0, ésspace (buf [p+ 4); F+4) { 

if (buf [p+ j] =?:’ V buf [p+ 7] =?1’) panic("Illegalcharacter,in,,itemname"); 
o, namebuf .str|j] = buf [p + J]; 


if (j =8A 7isspace (buf [p+ j])) panic("Item name, too, long"); 

00, Iname (last_itm < 2) = namebuf .lr.l, rname (lastitm < 2) = namebuf .lr.r; 

(Check for duplicate item name 16 ); 

last_itm-++; 

if (last_itm > maz_cols) panic("Too.many,items"); 

for (p += j +1; 0, isspace (buf [p]); p++) ; 

if (buf [p]= I?) { 
if (second 4 maz_cols) panic("Itemname,line,contains,,| twice"); 
second = last_itm, 
for (p++; 0, isspace (buf [p]); p++) ; 

i 

} 


This code is used in section 3. 


16. (Check for duplicate item name 16) = 
for (k = lastitm —1; k; k——) { 
if (0, Iname(k « 2) 4 namebuf.ir.l) continue; 
if (rname(k < 2) = namebuf.ir.r) break; 


if (k) panic("Duplicate,item,name" ); 


This code is used in section 15. 
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17. (Input the options 17) = 
while (1) { 
if (fgets (buf , bufsize, stdin)) break; 
line + ; 
if (0, buf [p = strlen(buf) — 1] 4 ?\n’) panic("Option line, too, long" ); 
for (p = 0; 0, isspace (buf [p]); p++) ; 
if (buf [p] =’ |? V abuf [p]) continue; /* bypass comment or blank line */ 
i = last_node; /*x remember the spacer at the left of this option */ 
for (pp = 0; buf[p]; ) { 
o, namebuf .Ir.l = namebuf .lr.r = 0; 
for (j =0; j < 8A (0, 7isspace (buf [p + j])) A buf [p+ j] 42:73 G++) 0, namebuf .str|j] = buf [p + J]; 
if (47) panic("Empty item name"); 
if (j =8A 7isspace (buf [p + j]) A buf [p+ 7] 4 ?:’) panic("Item mame, too,,long"); 
(Create a node for the item named in buf [p] 18); 
if (buf [p+ j] 4 ?:’) 0, nd[last_node].cilr = 0; 
else if (k > second) { 
if ((0, isspace (buf [p + 7 + 1])) V (0, misspace (buf [p + 7 + 2]))) 
panic("Color,must,be,a,single, character"); 
o, nd[last_node].clr = (unsigned char) buf [p+ 7 + 1); 
pt+=2; 
} else panic("Primaryitem must jbe,uncolored"); 
for (p += j +1; 0, isspace (buf [p]); p++) ; 


} 
if (spp) { 
if (vbose & show_warnings) fprintf (stderr, "Option, ignored,,(no.primary,jitems) :,"O"s", buf); 
while (last_node >i) { 
(Remove last_node from its item list 19); 


last_node —-; 
} else { 
o, nd[i].loc = last_node — 1; /* complete the previous spacer */ 
last_node ++; /* create the next spacer */ 
if (last_node = maz_nodes) panic("Too, many nodes"); 
options ++; 


o, nd[last_node].itm = i+ 1 — last_node; 


} (Initialize item 20); 
(Expand set 21); 
(Adjust nd 22); 


This code is used in section 3. 
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18. We temporarily use pos to recognize duplicate items in an option. 


(Create a node for the item named in buf [p] 18) = 
for (k = (last_itm — 1) « 2; kj k -=4) { 
if (0, Iname(k) A namebuf .lr.l) continue; 
if (rname(k) = namebuf .ir.r) break; 
} 
if (4k) panic("Unknown,item, name"); 
if (0, pos(k) > 7%) panic("Duplicate,item,mamein,,this,option"); 


last_node-++; 
if (last_node = maa_nodes) panic("Too many modes"); 
0, t = size(k); /* how many previous options have used this item? */ 


o, nd[last_node].itm = k >> 2, nd[last_node].loc = t, nd[last_node].linenr = line ; 
if ((k > 2) < second) pp =1; 
o, size(k) =t +1, pos(k) = last_node; 


This code is used in section 17. 


19. (Remove last_node from its item list 19) = 
0, k = nd[last_node].itm « 2; 
00, size(k)—-—, pos(k) =i-—1, 


This code is used in section 17. 


20. (Initialize item 20) = 
active = itemlength = last_itm — 1; 
for (k = 0,7 = primeztra; k < itemlength; k++) 
oo, item|k] = j,j += (k+2 < second ? primeztra : secondextra) + size((k +1) « 2); 
setlength = 7 — 4; /* a decent upper bound */ 
if (second = mazcols) osecond = active, second = j; 
else osecond = second — 1; 


This code is used in section 17. 
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21. Going from high to low, we now move the item names and sizes to their final positions (leaving room 


for the pointers into nb). 


(Expand set 21) = 
for (; k; k—) { 
0,j = item|k — 1); 
if (k = second) second = J; /* second is now an index into set */ 
00, size(j) = size(k < 2); 
if Gee j) =0Ak < osecond) baditem = k; 

0, ee j)=k-1, 

00, rname(j) = rnaiie (le < 2), Iname(j) = Iname(k « 2); 


} 


This code is used in section 17. 
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22. (Adjust nd 22) = 
for (k =1; k < last_node; k++) { 


if (0, nd[k].itm <0) continue; /* skip over a spacer */ 
0,j = ttem|nd|[k].itm — 1]; 
i= 7 + nd[k].loc; /* no mem charged because we just read nd|k].itm */ 


o, nd|k].ttm = j, nd[k].loc = t; 
o, set |i] = k; 
} 


This code is used in section 17. 


23. (Report an uncoverable item 23) = 
{ 
if (vbose & show_choices) { 
forintf (stderr, "Item"); 
print_item_name (item|baditem — 1], stderr); 
forintf (stderr, "\jhasyno,options!\n"); 


} 


This code is used in section 3. 


24. The “number of entries” includes spacers (because DLX2 includes spacers in its reports). If you want 
to know the sum of the option lengths, just subtract the number of options. 


(Report the successful completion of the input phase 24) = 
forintf (stderr, "("O"11ld options, ,"O"d+"O"d items, .,"O"d entries,successfullyread) \n", 
options , osecond, itemlength — osecond, last_node); 


This code is used in section 3. 


25. The item lengths after input are shown (on request). But there’s little use trying to show them after 
the process is done, since they are restored somewhat blindly. (Failures of the linked-list implementation in 
DLX2 could sometimes be detected by showing the final lengths; but that reasoning no longer applies.) 


(Report the item totals 25) = 
{ 

fprintf (stderr, "Item ,totals:"); 

for (k =0; k < itemlength; k++) { 
if (k = second) fprintf(stderr ,"|"); 
fprintf (stderr, "\\"O"d", size (item [k])); 

} 

fprintf (stderr, "\n"); 


} 


This code is used in section 3. 


26. (Randomize the item list 26) = 
for (k = active; k>1; ) { 
mems += 4,7 = gb_unif_rand(k); 
k--} 
00, 00,t = item|[j], item[j] = item|[k], item|k] = t; 
00, pos(t) = k, pos (item[j]) = j; 
} 


This code is used in section 3. 
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27. Proof logging. We will convert the XCC problem to a pseudo-boolean problem, so it can be verified 
by VeriPB. In the pseudo-boolean problem, all variables are booleans (0 or 1). We create one variable for 
each option (named Ojineny), and one variable for each item in each option. Primary items are named 
Olinenr-variablename, secondary items also include their color. 

Then we add constraints for each option/item pair, so that items are selected if and only if the option is 
selected. We add the constraint that for each primary item, the sum of all its corresponding ’item-in-option’ 
variables is 1, to make sure every primary item is covered. Finally, we add a constraint for each pair of 
incompatible options due to secondary items. 

We keep count of the number of constraints that we output during proof logging. Linear equalities count 
double, since they are converted to two linear inequalities by VeriPB. 

No mems are counted during proof logging. 

For example, let’s consider again the example from DLX2: 


| A simple example of color controls 
ABC|XY 

A B X:0 Y:0 

AC X:1 Y:1 

X:0 Y:1 

B X:1 

C Y:1 


The first option will be transformed in these 4 lines, to ensure the option and its variables are both either 
true or false: 


1 o38_A 1 ~03 = 1; 
1 0o38_B 1 ~03 = 1; 
1 038_X_0 1 ~03 = 1; 
1 o8_Y_0 1 ~03 = 1; 


To ensure item A will be covered, this line is added: 
1 o8_A 1 04_A = 1; 


Finally, to ensure item X isn’t given two different colors by the options of lines 3 and 4, we add this 
constraint: 
1 ~o3_X_0 1 ~04X_1 >= 1; 


(Output our problem to the OPB file for proof logging 27) = 
plconstraints = 0; 
(Output constraints for options and items 29); 
(Output constraints to ensure each primary item is covered 30); 
(Output constraints to ensure secondary items are assigned only one color 31); 
if (opb_file) fclose(opb_file); 


This code is used in section 3. 


28. (Subroutines 11) += 
void print_pLitemname(FILE «stream, int optnr,int k,int clr) 
{ 
namebuf .lr.l = Iname(k), namebuf .lr.r = rname(k); 
forintf (stream, "0o"O"a_"O" .8s", optnr, namebuf .str); 
if (k > second) { 
if (clr) fprintf (stream,"_"O"c", clr); 
else fprintf (stream,"_[]"); 
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29. (Output constraints for options and items 29) = 
for (1 =0; i < last_node; i++) { 
k = nd[i}.itm; 
if (k >0) { 
forintf (opb_file,"1u"); 
print_plitemname (opb-_file, nd|[2].linenr, k, nd{i].clr); 


forintf (opb_file, "41470" O"dy=u13;\n", nd[i].linenr, namebuf .str, optnr); 
plconstraints += 2; 


} 
} 


This code is used in section 27. 


30. 


{ 


(Output constraints to ensure each primary item is covered 30) = 


int aw, ss, opt; 
for (k =0; k < active; k++) { 
ii = item|k]; 
if (ii > second) continue; 
for (s = ti, ss = s+ size(w#); s < ss; st+) { 
opt = set|s]; 
printf (opb_file,"14."); 
print_pLitemname ( opb_file, nd[opt].linenr , item |[k], 0); 
{printf (opb_file, ""); 
} 
forint (opb_file, "=,1;\n"); 
pl.constraints += 2; 
} 


} 


This code is used in section 27. 
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31. For the secondary items, we add a condition in the OPB file for each pair of incompatible options. 


(Output constraints to ensure secondary items are assigned only one color 31) = 


{ 
int w, s1, s2, ss, optl, opt2, opta; 
for (k =0; k < active; k++) { 
ii = item|k]; 
if (ii < second) continue; 
for (si = ti,ss = s1 + size(ti); s1 < ss; s1++) { 
opt1 = set|s1]; 
for (s2 =s1 +1; s2 < ss; s2++) { 
opt2 = set[s2]; 
if (nd[opt1].clr 4 nd[opt2].clr V nd[opt1].clr =0V nd[opt2].clr =0) { 
fprintf (opb_file, "147"); 
print_pLitemname ( opb_file, nd[opt1 ].linenr, item|[k], nd [opt |.clr); 
forintf (opb_file, "41u~"); 
print_pLitemname ( opb_file, nd[opt2].linenr, item[k], nd|[opt2|.clr ); 
forintf (opb_file, "4>=1;\n"); 
plconstraints + ; 


This code is used in section 27. 


32. We also need to write the necessary steps to the OPB file, so that VeriPB can verify our proof. This 
uses the routines below. 
(Output header of the PBP file for proof logging 32) = 

forint (pbp_file, "pseudo-Boolean, proof,,version,1.0\n"); 

forintf (pbp_file, "£,"O"d\n", pl_constraints ); 


This code is used in section 3. 


33. (Write a solution for proof logging 33) = 
sprintf (pbp_file, "vu"); 
int opt; 
for (k =0; k < level; k++) { 
for (opt = choice[k]; ndlopt].itm > 0; opt——) ; 
forintf (pbp_file, "o"O"d.", nd[+ opt].linenr); 
for (opt; nd[opt].itm > 0; opt++) { 
print_plitemname (pbp-file, nd[opt].linenr, nd[opt].itm, nd[opt].clr); 
Sprintf (pbp-file, "u"); 


} 
printf (pbp_file,"\n"); 
pl_constraints ++; 


This code is used in section 48. 
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34. When backtracking, we need to add a RUP constraint. 
(Write RUP constraint on backtrack for proof logging 34) = 


{ 
int opt; 
if (level > 0) { 
Sprintf (pbp-file, "rup.s" ); 
for (k =0; k < level; k++) { 
fprintf (pbp_file, "147 0"O"d,", nd[ choice [k]].linenr); 
} 
forintf (pbp_file, ">=,1;\n"); 
pl_constraints ++; 
} 
} 


This code is used in section 38. 


35. When hiding failes, we also need to add a RUP constraint. 
(Write RUP constraint on failed hide for proof logging 35) = 


{ 
int opt; 
if (level > 0) { 
Sprintf (pbp-file, "rup." ); 
for (k =0; k < level; k++) { 
fprintf (pbp_file, "147 0"O"d.", nd[ choice [k]].linenr); 


fprintf (pbp_file, "1,~0"O"d>=y1 ;\n", nd[cur_node].linenr); 
pl_constraints ++; 
} 


} 


This code is used in section 43. 
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36. When we finish our program, we need to notify our proof logging that we should have derived a 


contradiction. 


( Write contradiction line for proof logging 36) = 
printf (pbp_file, "cu-1\n"); 


This code is used in section 3. 


37. For better performance of the verifier, it is useful to remove derived constraints from deeper levels of 
our backtrack tree, once they become reduntant for the further verification. For each level of our backtrack 
tree, we log when we enter this level. When we backtrack, we wipe all constraints derived on the lower levels, 


since these have become reduntant. 


(Subroutines 11) += 
void write_proof_log_level (int level) 


fprintf (pbp_file, "#4"O"d\n", level); 
void wipe_proof_log_level (int level) 


forintf (pbp_file, "w4"O"d\n", level); 
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38. The dancing. Our strategy for generating all exact covers will be to repeatedly choose an item that 
appears to be hardest to cover, namely an item whose set is currently smallest, among all items that still 
need to be covered. And we explore all possibilities via depth-first search. 

The neat part of this algorithm is the way the sets are maintained. Depth-first search means last-in-first- 
out maintenance of data structures; and the sparse-set representations make it particularly easy to undo 
what we’ve done at deeper levels. 

The basic operation is “covering an item.” That means removing it from the set of items needing to be 
covered, and “hiding” its options: removing them from the sets of the other items they contain. 


(Solve the problem 38) = 
{ 
level = 0; 
forward: nodes ++; 
write_proof_log_level (level ); 
if (vbose & show-_profile) profile |level]++; 
if (sanity_checking) sanity(); 
(Maybe do a forced move 47); 
(Do special things if enough mems have accumulated 40 ); 
(Set best_itm to the best item for branching 46 ); 
if (forced) { 
o, best_itm = force|—— forced]; 
(Do a forced move 51); 


if (¢ = maz_nodes) (Visit a solution and goto backup 48); 
(Swap best_itm out of the active list 41); 
oactive = active, hide (best_itm, 0,0); /* hide its options */ 
cur_choice = best_itm; 
(Save the currently active sizes 49); 
advance: oo, cur.node = choice [level] = set[cur_choice]; 
tryit: if ((vbose & show_choices) A level < show_choices_max) { 
forintf (stderr, "L"O"d:", level); 
print_option(cur_node, stderr); 


} 
(Swap out all other items of cur_node 42); 
(Hide the other options of those items, or goto abort 43); 
if (++level > maal) { 
if (level > maz_level) { 
forintf (stderr, "Too,manyjlevels!\n"); 
exit (—4); 


maal = level; 

} 

goto forward; 
backup: 

if (level > 0) write_proof_log_level (level — 1); 

(Write RUP constraint on backtrack for proof logging 34); 

if (level > 0) wipe_proof_log_level (level ); 

if (level =0) goto done; 

level —; 

00, cur_node = choice [level], best_itm = nd|cur_node].itm, cur_choice = nd[cur_node].loc; 
abort: if (0, cur_choice +1 > best_itm + size(best_itm)) goto backup; 

(Restore the currently active sizes 50); 

cur_choice++; goto advance; 
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} 


This code is used in section 3. 


39. We save the sizes of active items on savestack, whose entries have two fields / and r, for an item and 
its size. This stack makes it easy to undo all deletions, by simply restoring the former sizes. 
(Global variables 4) += 

int level; /* number of choices in current partial solution */ 

int choice |maz_level]; /* the node chosen on each level */ 

int saved|maz_level + 1]; /* size of savestack on each level */ 

ullng profile |maz_level]; /* number of search tree nodes on each level */ 

twoints savestack |savesize]; 


int saveptr; /* current size of savestack */ 


40. (Do special things if enough mems have accumulated 40) = 
if (delta A (mems > thresh)) { 
thresh += delta; 
if (vbose & show_fulLstate) print_state(); 
else print_progress(); 
} 
if (mems > timeout) { 
fprintf (stderr, "TIMEOUT! \n"); goto done; 


This code is used in section 38. 


41. (Swap best_itm out of the active list 41) = 
p = active — 1, active = p; 
0, pp = pos(best_itm); 
0, cc = ttem|[p]; 


00, item |p] = best_itm, item[pp] = cc; 
00, pos(cc) = pp, pos(best_itm) = p; 
updates ++; 


This code is used in sections 38 and 51. 
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42. Note that a colored secondary item might have already been purified, in which case it has already been 
swapped out. We don’t want to tamper with any of the inactive items. 
(Swap out all other items of cur_node 42) = 
p = oactive = active; 
for (q = cur_node +1; q # cur_node; ) { 
0,c = nd|q].itm; 
if (<0) ¢q+=c¢ 


else { 
0, pp = pos(c); 
if (pp <p) { 


0, cc = item|——p]; 

00, item |p] = c, item|[pp] = cc; 
00, pos(cc) = pp, pos(c) = p; 
updates ++; 


q+; 
} 
} 
active = p; 


This code is used in section 38. 


43. A secondary item was purified at lower levels if and only if its position is > oactive. 


(Hide the other options of those items, or goto abort 43) = 
for (q = cur_node +1; q #4 cur_node; ) { 

0, cc = nd[q].itm; 

if (cc <0) q+= cc; 

else { 

if (cc < second) { 
if (hide(cc,0,1) =0) { 

forced = 0; 
(Write RUP constraint on failed hide for proof logging 35 ); 
goto abort; 


} else { /* do nothing if cc already purified «/ 
0, pp = pos(cc); 
if (pp < oactive A (0, hide (cc, nd|q}.clr,1) =0)) { 
forced = 0; 
(Write RUP constraint on failed hide for proof logging 35 ); 
goto abort; 


This code is used in section 38. 
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44. The hide routine hides all of the incompatible options remaining in the set of a given item. 

If check is nonzero, further checking is done: (1) hide returns zero if its actions would cause a primary 
item to be uncoverable. (2) An active primary item that is now coverable in only one way is placed on the 
force stack. 

If the color parameter is zero, all options are incompatible. Otherwise, however, the given item is 
secondary, and we retain options for which that item has a color match. 

When an option is hidden, it leaves all sets except the set of that given item. And the given item is 
inactive. Thus a node is never removed from a set twice. 


(Subroutines 11) += 
int hide(int c,int color,int check) 
{ 
register int cc, s, rr, ss, nn, tt, uu, vv, nnp; 
for (0, rr =c,s =c+ size(c); rr <8; rr++) { 
o, tt = set [rr]; 
if (>color V (0, nd[tt].clr 4 color)) (Remove option tt from the other sets it’s in 45); 


} 


return 1; 


} 


45. (Remove option tt from the other sets it’s in 45) = 
{ 
for (nn = tt +1; nn Att; ) { 
0, uu = nd[nn].itm, vv = nd[nn].loc; 
if (uu <0) { nn += uu; continue; } 
if (0, pos(uu) < oactive) { 
0, 88 = size(uu) — 1; 
if (ss <1A check A uu < second A pos(uu) < active) { 
if (ss =0) { 
if ((ubose & show_choices ) A level < show_choices.max) { 
forintf (stderr, "\ycan’ tycover"); 
print_item_name(uu, stderr ); 
forintf (stderr, "\n"); 
} 
return 0; 
} else 0, force[forced++] = uu; 
} 
0, nnp = set[uu + ss]; 
0, size(uu) = ss; 
oo, set[uu + ss] = nn, set[vu] = nnp; 
oo, nd[nn].loc = uu + ss,nd[nnp].loc = vv; 
updates ++; 
} 
nn+; 
} 
} 


This code is used in section 44. 
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46. The “best item” is considered to be an item that minimizes the number of remaining choices. All 
candidates of size 1, if any, are put on the force stack. If there are several candidates of size > 1, we choose 
the leftmost. 

Notice that a secondary item is active if and only if it has not been purified (that is, if and only if it hasn’t 
yet appeared in a chosen option). 

(This program explores the search space in a slightly different order from DLX2, because the ordering of 
items in the active list is no longer fixed. But ties are broken in the same way when s > 1.) 


(Set best_itm to the best item for branching 46) = 
t = max_nodes , tmems = mems; 
if ((vbose & show_details) A level < show_choices_max / level > maal — show_choices_gap ) 
fprintf (stderr, "Level"O"d:", level); 
for (k =0; k < active; k++) 
if (0, item[k] < second) { 
0, 8 = size(item|k]); 
if ((vbose & show_details) \ level < show_choices_max / level > maal — show_choices_gap) { 
print_item_name (item|[k], stderr); 
fprintf (stderr, "("O"d)", 8); 


} 
if (s <1) { 
if (s =0) fprintf(stderr, "I’m confused. \n"); /* hide missed this */ 
else o, force|forced ++] = item|k]; 
} else if (s <t) { 
if (s < t) best_itm = item|k],t = s; 
else if (item[k] < best_itm) best_itm = item|k]; /* suggested by P. Weigel */ 


} 


if ((vbose & show_details) \ level < show choices.max / level > maal — show-choices.gap) { 
if (forced) fprintf(stderr, "found,"O"d forced\n", forced); 
else if (¢ = maz_nodes) fprintf(stderr, "jsolution\n"); 
else { 
fprintf (stderr, "\joranching,,on"); 
print_item_name (best_itm, stderr ); 
forintf (stderr, "("O"d)\n",t); 


} 


if (t > mardeg At < maz_nodes \ forced) maxdeg = ¢; 
if (shape_file) { 
if (¢ = maz_nodes) fprintf (shape_file, "sol\n"); 
else { 
forintf (shape_file,""O"d", t); 
print_item_name (best_itm, shape_file); 
forintf (shape_file, "\n"); 


fflush(shape_file); 
} 


cmems += mems — tmems; 


This code is used in section 38. 
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47. (Maybe do a forced move 47) = 
while (forced) { 
o, best_itm = force|—— forced]; 
if (0, pos(best_itm) < active) (Do a forced move 51); 


} 


This code is used in section 38. 


48. (Visit a solution and goto backup 48) = 
4 

count +; 

(Write a solution for proof logging 33); 

if (spacing A (count mod spacing =0)) { 
printf (""O"11d:\n", count); 
for (k =0; k < level; k++) print_option (choice |k], stdout); 
fflush (stdout); 


if (count > maxcount) goto done; 
goto backup; 


} 


This code is used in section 38. 


49. (Save the currently active sizes 49) = 
if (saveptr + active > maxsaveptr) { 
if (saveptr + active > savesize) { 
forintf (stderr, "Stack overflow, (savesize="O"d) !\n", savesize); 


exit (—5); 
} 
maxsaveptr = saveptr + active; 
} 
for (p = 0; p < active; pt+) 


000, savestack [saveptr + p].l = item|[p], savestack |saveptr + p|.r = size (item [p}); 
o, saved [level + 1] = saveptr = saveptr + active; 


This code is used in section 38. 


50. (Restore the currently active sizes 50) = 
0, saveptr = saved |level + 1]; 
0, active = saveptr — saved [level]; 
for (p = —active; p <0; pt+) 00, size(savestack |saveptr + p].l) = savestack|saveptr + p].r; 


This code is used in section 38. 
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51. <A forced move occurs when best_itm has only one remaining option. In this case we can streamline 


the computation, because there’s no need to save the current active sizes. (They won’t be looked at.) 
(Do a forced move 51) = 


{ 
if ((ubose & show_choices) A level < show_choices.maz) fprintf(stderr," (forcing) \n"); 
(Swap best_itm out of the active list 41); 
oactive = active, hide (best_itm, 0,0); /* hide its options */ 
cur_choice = best_itm; 
o, saved [level + 1] = saveptr; /* nothing placed on savestack */ 
goto advance; 


} 


This code is used in sections 38 and 47. 
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52. (Subroutines 11) += 
void print_savestack (int start,int stop) 
{ 
register k; 
for (k = start; k < stop; k++) { 
print_item_name (savestack [k].1, stderr); 
forintf (stderr, "("O"d) ,u"O"d\n", savestack [k].1, savestack [k].r); 


} 


53. (Subroutines 11) += 
void print_state (void) 
{ 
register int 1; 
forintf (stderr, "Currentystate,(level,"O"d) :\n", level); 
for (J =0; | < level; I++) { 
print_option (choice [l], stderr ); 
if (1 > show_levels_maz) { 
forintf (stderr, "\...\n"); 
break; 
} 
i 


forintf (stderr, "\y"O"1ld solutions, ,"O"1ld.mems , ,and,max,jlevel,"O"dsoyfar.\n", count, 


mems, mazl); 
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54. During a long run, it’s helpful to have some way to measure progress. The following routine prints a 
string that indicates roughly where we are in the search tree. The string consists of character pairs, separated 
by blanks, where each character pair represents a branch of the search tree. When a node has d descendants 
and we are working on the kth, the two characters respectively represent k and d in a simple code; namely, 
the values 0, 1, ..., 61 are denoted by 


All values greater than 61 are shown as ‘*’. Notice that as computation proceeds, this string will increase 
lexicographically. 

Following that string, a fractional estimate of total progress is computed, based on the naive assumption 
that the search tree has a uniform branching structure. If the tree consists of a single node, this estimate 
is .5; otherwise, if the first choice is ‘k of d’, the estimate is (k —1)/d plus 1/d times the recursively evaluated 
estimate for the kth subtree. (This estimate might obviously be very misleading, in some cases, but at least 
it tends to grow monotonically.) 


(Subroutines 11) += 
void print_progress (void) 
{ 
register int 1, k, d, c, p, ds =0; 
register double f/f, jd; 
fprintf (stderr, "\after,"O"1lld mems:"O"1lldsols,", mems, count); 
for (f = 0.0, fd = 1.0,1 =0; I < level; 1++) { 
c= nd{choice[l]].itm,d = size(c),k = nd{choice[l]].loc —c +1, 
fd x=d, f += (k —1)/fd; /* choice | is k of d */ 
if (l < show_levels_maz ) 
forintf (stderr, "W"O"c"O"c",k <10??0°+k:k < 36? ’a’+k—-10:k < 62? °A? +k—36: ?*?, 
d<10??0? +d:d<36?’a’+d—10:d< 62? 7A’? +d—36: ’*?’); 
else if (ads) ds = 1, fprintf (stderr,"..."); 


} 
forintf (stderr, "\"O" .5£\n", f + 0.5/fd); 
} 


55. (Print the profile 55) = 


forintf (stderr, "Profile:\n"); 
for (level = 0; level < mazal; level++) fprintf(stderr ,""O"3d:"O"1ld\n", level, profile [level]); 
} 


This code is used in section 3. 
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cur_node: 3, 35, 38, 42, 43. 


ae 4, 5, 40. ee : ef 11, 15, 16, 17, 18, 28, 29 
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plconstraints: 4, 27, 29, 30, 31, 32, 33, 34, 35. stderr: 2, 3, 4, 5, 12, 18, 14, 15, 17, 23, 24, 25, 
pos: 8, 13, 14, 18, 19, 21, 26, 41, 42, 43, 45, 47. 38, 40, 45, 46, 49, 51, 52, 53, 54, 55. 
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print_pLitemname: 28, 29, 30, 31, 33. strepy: 5. 
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print_state: 40, 53. strlen: 15, 17. 
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prow: 12, 13. t: 3. 
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rr: 44. uint: 3. 
3: 3, 44. ullng: 3, 4, 39. 
sanity: 14, 38. updates: 3, 4, 41, 42, 45. 
sanity_checking: 3, 14, 38. uu: 44, 45. 
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secondextra: 8, 20. 

set: 8,9, 10, 18, 14, 15, 21, 22, 30, 31, 38, 44, 45. 

setlength: 3, 9, 13, 20. 

shape_file: 4, 5, 6, 46. 

shape.name: 4, 5. 

show_basics: 3, 4. 

show-choices: 4, 23, 38, 45, 51. 

show_choices_gap: 4, 5, 46. 

show_choices.max: 4, 5, 38, 45, 46, 51. 

show_details: 4, 46. 

show-_full_state: 4, 40. 

show_levels.maz: 4, 5, 53, 54. 

show_maxz_deg: 3, 4. 

show-profile: 3, 4, 38. 

show_tots: 3, 4. 

show.warnings: 4, 17. 

size: 8, 12, 13, 14, 18, 19, 20, 21, 25, 30, 31, 
38, 44, 45, 46, 49, 50, 54. 

spacing: 4, 5, 48. 

ss: 30, 31, 44, 45. 

sscanf: 5. 

start: 52. 
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(Adjust nd 22) Used in section 17. 

(Check for duplicate item name 16) Used in section 15. 
(Close the files 6) Used in section 3. 

(Create a node for the item named in buf [p] 18) Used in section 17. 
(Do a forced move 51) Used in sections 38 and 47. 
(Do special things if enough mems have accumulated 40) Used in section 38. 
(Expand set 21) Used in section 17. 

(Global variables 4, 9, 39 ) Used in section 3. 
(Hide the other options of those items, or goto abort 43) Used in section 38. 
(Initialize item 20) Used in section 17. 

(Input the item names 15) Used in section 3. 
(Input the options 17) Used in section 3. 
(Maybe do a forced move 47 ) Used in section 38. 
(Output constraints for options and items 29) Used in section 27. 
(Output constraints to ensure each primary item is covered 30) Used in section 27. 
(Output constraints to ensure secondary items are assigned only one color 31) Used in section 27. 
(Output header of the PBP file for proof logging 32) Used in section 3. 
(Output our problem to the OPB file for proof logging 27) Used in section 3. 
(Print the profile 55) Used in section 3. 
(Process the command line 5) Used in section 3. 
(Randomize the item list 26) Used in section 3. 

(Remove option ét from the other sets it’s in 45) Used in section 44. 

(Remove last_node from its item list 19) Used in section 17. 

(Report an uncoverable item 23) Used in section 3. 

(Report the item totals 25) Used in section 3. 

(Report the successful completion of the input phase 24) Used in section 3. 
(Restore the currently active sizes 50) Used in section 38. 

(Save the currently active sizes 49) Used in section 38. 

(Set best_itm to the best item for branching 46) Used in section 38. 

(Solve the problem 38) Used in section 3. 

(Subroutines 11, 12, 13, 14, 28, 37, 44, 52, 53,54) Used in section 3. 

(Swap out all other items of cur_node 42) Used in section 38. 

(Swap best_itm out of the active list 41) Used in sections 38 and 51. 

(Type definitions 8,10) Used in section 3. 

( Visit a solution and goto backup 48) Used in section 38. 

(Write RUP constraint on backtrack for proof logging 34) Used in section 38. 
(Write RUP constraint on failed hide for proof logging 35) Used in section 43. 
( Write a solution for proof logging 33) Used in section 48. 

( Write contradiction line for proof logging 36) Used in section 3. 
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